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Inspired by recent experiments where electron transport was measured across graphene nanorib¬ 
bons (GNR) suspended between a metal surface and the tip of a scanning tunneling microscope 
[Koch et al, Nat. Nanotechnol. 7, 713 (2012)], we present detailed first-principles simulations of 
inelastic electron tunneling spectroscopy (lETS) of long pristine and defected armchair and zigzag 
nanoribbons under a range of charge carrier conditions. For the armchair ribbons we find two robust 
lETS signals around 169 and 196 mV corresponding to the D- and G-modes of Raman spectroscopy 
as well as additional fingerprints due to various types of defects in the edge passivation. For the 
zigzag ribbons we show that the spin state strongly influences the spectrum and thus propose lETS 
as an indirect proof of spin polarization. 

PACS numbers: 81.05.ue, 73.63.-b, 72.10.Di 


I. INTRODUCTION 

Graphene as the basis of a new generation of 
electronically has been the center of much attention in 
the last years, and devices based on nanostructured 
graphene have been put forward. The most generic form 
of nanostructured graphene is graphene nanoribbons 
(GNR)^and other structures, such as graphene anti-dot 
lattices^, can be viewed as networks of them. GNRs 
are potential candidates for molecular wires with tailored 
conductance properties. For graphene-based nanostruc¬ 
tures the edges and their passivation, as well as defects 
inside the structure, can play crucial roles for the trans¬ 
port properties.iy However, characterization of edge passi¬ 
vation or structural/chemical defects is challenging espe¬ 
cially after device fabrication. Raman spectroscop}!^ can 
give information about defects on large areas of the sam¬ 
ple, while tip-enhanced Raman spectroscopy (TERS)|y 
have been used in combination with STM on GNRs. 
However, Raman studies involve averages over larger ar¬ 
eas (> 10 nm), and does not yield information about 
the impact of vibrations on transport. In that aspect 
inelastic electron tunneling spectroscopy (lETS) serves 
as a way of performing non-destructive characterization 
yielding vibrational/phonon fingerprints of a range of de¬ 
fects. In order to interpret lETS experiments, theoretical 
modeling of the inelastic signals in the electronic current 
due to electron-phonon (e-ph) scattering is needed. 

GNRs have been fabricated using different strate¬ 
gies inc ludin g lithographic techniques ,1^ chemical 
synthesisepitaxial growtli^^, and longitudinal 
unzipping of carbon nanotubes.^ Furthermore, several 
groups have succeeded in atomically preci se bo ttom-np 
fabrication of armchair GNRs (AGNR)P3III1^ chiral 
GNRs,^ and AGNR hetero-junction^^ grown on metal 
surfaces. Experimentally, the vibrational properties 
have been investigated by Raman spectroscopy and 


the electronic structure has been mapped out by 
STM, angle-resolved (two-photon) photo-emissio n and 
high-resolution electron energy loss spectroscopyj ^ l ^^ l ^ y 
Signatures of phonon excitation were observed by STM 
in the differential conductance spectroscopy performed 
at the zigzag termini state of AGNRs adsorbed on 
Au(lll), and these signatures were shown to be sen¬ 
sitive to modifications in the local atomic geometry.l^ 
AGNRs have also been lifted up from the weakly 
bonding Au(lll) surface with the tip of a STM enabling 
measurements of the vol tag e-dependent conductance in 
suspended configurations.^ 

From the theoretical side density-functional the¬ 
ory (DFT) has been used to investigate the stabil¬ 
ity of structural and chemical reconstructions of GNR 
edgesp^Hm tog ether with the transport and band- 
gap engineeringxhe vibrational properties and 
phonon band structure hav e bee n calculated with empir¬ 
ical potentially and DFT.^^^^ In addition, there have 
been theoretical predict ion^^y^ of the Ram an spectrum, 
in good agreement with experimentsFor a finite 
AGNR the role of zigzag termini states have been studied 

DFT to the many-body Hub- 

Inspired by the recent lifting experiments by Koch 
et al. ^ we here investigate theoretically the signals of e- 
ph scattering in the conductance of long GNRs between 
metal electrodes. Our aim is two-fold. First, we want to 
address the role phonon scattering in the transport char¬ 
acteristics of pristine GNRs. Second, we wish to compute 
detailed lETS for different GNRs under varying charge 
carrier concentrations and explore how different types of 
realistic defects may modify the lETS and thus possibly 
be directly probed in transport measurements. We focus 
on the two most generic edge types, namely armchair 
(AGNR) and zigzag (ZGNR), and pay attention to the 
effects of spin polarization in the latter case. In actual 


theoretically, comparing 
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EIG. 1: (Color online) (a) Computational setup for a pristine AGNR showing electrode, device and dynamical regions, (b) 
Electronic band structure {kx is in units of inverse unit cell length). The different bands are colored according to symmetry 
of the electronic states. Red: symmetric, corresponding to Eig. [^a-b). Blue: anti-symmetric, corresponding to Eig. Hc-d). 
(c) Electronic transmission for varying electrode broadening describing the coupling to the metal contacts, 77 = 0, 0.1,1 eV, see 
text, (d) Electronic DOS projected onto the dynamical region. Panels (e)-(h) show the similar entities for the pristine ZGNR 
case. 
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EIG. 2: (Color online) (a)-(d) Electron transmission eigen- 
channels for the clean AGNR for the valence bands at — 
Ep = —1 eV and for the conduction bands Sii E — Ep = 1 
eV. (e)-(h) Electron transmission eigenchannels for the clean 
ZGNR in the valence bands at E — Ep = —0.4 eV and in the 
conduction bands at E — Ep = 0.4 eV for one spin component. 
The eigenchannels for the other spin component are simply 
mirror images around the middle of the ZGNR (not shown). 
The red/blue (pink/gray) isosurfaces represent the real (imag¬ 
inary) part and sign of the scattering state wave function. Eor 
all eigenchannel calculations the electrode broadening was set 
to zero (?7 = 0 eV). 


experiments the substrate or an applied gate potential 
control the Fermi level Ep in the ribbons. To address 
this variability we scan Ep using a numerically effective 
scheme enabling fast calculations of the lETS.^ We find 
that the AGNR generally display two robust lETS sig¬ 
nals around 169 and 196 mV corresponding to the D- 
and G-modes of Raman spectroscopy and that a dehy¬ 
drogenated dimer at the edge should further leave a clear 
defect signal at around 245 mV. Eor the ZGNR we find 
that the spin polarization breaks the mirror symmetry 
around the middle of the ribbon resulting in lETS sig¬ 
nals from a range of modes around the D- and G-mode 
energies. Eor both AGNR and ZNGR defects which break 
the planar symmetry of ribbons allows for contributions 
to the lETS from out-of-plane phonon modes. 


The paper is organized as follows. Eirst we discuss our 
atomistic model setup for the density functional and elec¬ 
tron transport calculations, and outline the approach for 
the lETS simulations. In Sec. Ill we present our results 
for pristine AGNR and ZGNR and relate their transport 
properties and lETS to the band structures. In Sec. IV 
we turn to the defected systems by considering realistic 
possibilities of defects in the edge passivation, backbone 
bonding motifs, and presence of adatoms. Einally, a sum¬ 
mary and our conclusions are presented in Sec. V. 
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II. METHODS 


We calculate the electronic and vibrational 
structure from D FT u sing the academic codes 
Siesta/TranSiestaP^I^ We employ the generalized 
gradient approximation (GGA) for the exchange- 
correlation functionalp^ a single-zeta polarized (SZP) 
basis set for the carbon and hydrogen atoms, and use 
a cut-off energy of 400-500 Ry for the real-space grid. 
These choices, balancing accuracy and computational 
cost, provide a good description to investigate trends and 
general behavior of the substantial number of systems 
considered in this work. 

The vibrational degrees of freedom, calculated by di- 
agonalization of the dynamical matrix, and the e-ph 
couplings are extracted from finite differences as imple¬ 
mented in the Inelastica code.^^^®^ The armchair and 
zigzag GNRs considered here are shown in Fig. We 
adopt the usual two-probe setup with the device region 
(D) coupled to left (L) and right (R) electrodes with 
all electronic matrix elements expressed in a local basis 
set. The primitive unit cell of the AGNR (ZGNR) con¬ 
sists of 18 (10) atoms and in our calculations this unit 
cell is repeated 10 (18) times in the transport direction 
to form the scattering regions illustrated in Fig. ^a,e). 
The electrode couplings are included on the two 

first/last unit cells before folding onto D. In our treat¬ 
ment a subset of atoms in D is allowed to vibrate. We 
fix this dynamical region, restricted by the condition that 
the e-ph couplings are fully included inside D, to the 4 
and 6 central unit-cells for the AGNR and ZGNR, respec¬ 
tively. The corresponding e-ph couplings used to calcu¬ 
late the inelastic electron transport are thus expressed in 
the center 6 unit-cells for the AGNR and 8 unit-cells for 
the ZGNR. The convergence of our results with the size 
of the dynamical region is addressed below. 

We generally consider nanoribbons that are suspended 
between two metallic leads. In the case of the lifting 
exp eriments,!^ these would correspond to the metal sam¬ 
ple surface and the STM tip. We wish here to focus on 
the action inside the GNRs and put aside the possible 
complications due to the detailed electronic structure of 
the metals, and the metal-GNR interface in particular. 
To this end we introduce a simple model of the metal 
electrodes without substantial electronic features: we use 
semi-infinite GNRs with highly broadened states (effec¬ 
tively smearing out energy gaps). In practice this is done 
by adding a finite numerical imaginary part r] to the en¬ 
ergy argument in the electrode recursion calculation.^^ 
This scheme ensures that the phonon effects originate 
from the GNRs themselves and not from details of the 
metal-GNR interface, which is generally unknown in the 
STM experiments. The electronic band structures for the 
infinite ribbons, along with the transmission and den¬ 
sity of states (DOS) are shown for r] = 0,0.1,1 eV in 
Fig. [3b,c,d) and Fig. [3f,g,h) for AGNR and ZGNR, 
respectively. We note that the broadened transmission 
spectrum [Fig. [^d)] is quite consistent with the exper¬ 


imentally reported differential conductance curves for 
AGNR.^ The electronic states involved in the trans¬ 
port are shown in Fig. in terms of the transmission 
eigenchannelJ^ in the valence and conduction bands of 
the AGNR and ZGNR. Their spatial symmetry play a 
significant role for the selection rules involved in the in¬ 
elastic scattering as discussed later. 

In principle, the electronic structure should be evalu¬ 
ated at finite bias. However, without a detailed model 
of the connection to the metal electrodes (where an im¬ 
portant part of the voltage drop will take place) and for 
sufficiently long systems (in which the electric field will be 
small), it is reasonable to use the zero-voltage electronic 
structure and to simply assume a symmetric voltage drop 
over the two identical, idealized device-electrode inter¬ 
faces. More specifically, in the following we consider that 
the chemical potentials of the electrodes move according 
to the applied bias voltage V and possibly an applied gate 
voltage Vg (mimicking actual doping or electrostatic gat¬ 
ing that modify the charge carrier concentration in an ex¬ 
perimental setup) according to = F’f +Vb- 


A. Computational scheme for lETS 


For a device strongly coupled to the electrodes, a cou¬ 
pling between the electron current I{V) and a phonon 
mode A ideally shows up at zero temperature as a step 
discontinuity in the differential conductance when the in¬ 
elastic phonon emission process becomes energetically al¬ 
lowed, that is, when the chemical potential difference ex¬ 
ceeds the quantum of vibrational energy, |/iL ~ = 

hwx. Thus, around the emission threshold the electronic 
states involved in the scattering process are those at 
and jj^R. The lETS signal, conventionally expressed as 
the ratio between the second and first derivatives of the 
current with respect to the voltage. 


lETS = 


dim 

dvm ’ 


( 1 ) 


is calculated by considering the e-ph coupling as the per¬ 
turbation on the current, evaluated using the nonequilib¬ 
rium Green’s functions (NEGE). In the so-called lowest 
order expansion (LOE) the inelastic part of the differen¬ 
tial conductance can be written 


dvHV) =7A dvTy^iV, T) (2) 

where summation over the vibration index A is assumed. 

and are the “universal” (system-independent) 

functions that depend on the applied bias E, phonon 
energy hwx and the temperature T. Assuming the elec¬ 
tronic and phononic distribution functions are given by 
the Eermi-Dirac and Bose-Einstein distributions, respec- 
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lively, their analytical expressions can be written as: 
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( 3 ) 


( 4 ) 


where Go = 2e^/h is the conductance quantum, f{e) is 
the Fermi-Dirac function, = 5 ' + shwx, and H denotes 
the Hilbert transform. 

The signal amplitudes 7 a and hzx of the symmetric and 
antisymmetric signals in the differential conductance are 
even and odd in bias, respectively. For a symmetric struc¬ 
ture the asymmetric signal vanishes in the wide-band ap¬ 
proximation (LOE-WBA^^. However, this is not guar¬ 
anteed in the more general treatment employed here,l^ 
where the energy dependence of the electronic structure 
is explicitly taken into account. The amplitudes 7 a and 
Kx are expressed in terms of electronic structure quanti¬ 
ties and e-ph couplings ,1^ 
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FIG. 3: (Color online) Convergence of the intrinsic lETS for 
pristine (a) AGNR and (b) ZGNR as a function of the size 
of the dynamical region (stated in the legends). The results 
are normalized with respect to the number of vibrating unit 
cells, i.e., we show the lETS amplitude per H 4 C 14 segment 
for AGNR and per H 2 C 8 segment for ZGNR. No gate voltage 
is applied (Vg = 0.0 V). 
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III. PRISTINE GRAPHENE NANORIBBONS 


7a =Tr[MAAL(/iL)MAAi?(/ii?)] + Im^A, (5) 

Kx =2ReBx, (6) 

where Bx is defined as 

Bx =Tr[MAA/,(/i^)r^(/i^)G’^(/i^)MAA/,(/r/,) 

- MxG^{/j.R)TL{fJ^R)AR{/j.R)MxAL{fiL)]^ (7) 

In the above, Ma denotes the e-ph coupling matrix 
for mode A, the retarded/advanced unperturbed 

Green’s functions, and A^^ = the spectral den¬ 

sity matrices for left/right moving states with the time- 
reversed version The purely electronic 

quantities are thus being evaluated at the chemical po¬ 
tentials of the left/right electrodes corresponding to the 
excitation threshold for each vibration. We compute Ma 
with the finite-difference scheme of Inelastica taking 
the vacuum energy as a common reference(in absence of 
real metal leads to pin the Eermi energy) 

In the localized atomic basis set of Siesta all the above 
quantities are matrices defined in the electronic space 
corresponding to region D. The second derivatives of 
the universal functions in Eqs. (§-0 are sharply peaked 
around the phonon threshold. Eor this reason the coeffi¬ 
cients 7a and hzx can be considered volt age-independent 
with their values computed exactly at the threshold. Due 
to the computational efficiency of the LOE scheme de¬ 
scribed above we are able to evaluate the lETS on a fine 
grid of gate voltages Vg spanning a large range of rele¬ 
vant values between valence and conduction bands of the 
GNRs. 


Now we first turn to the lETS results of the two pris¬ 
tine (clean) ribbons, and in the following section to the 
impact of selected defects in the lETS. As our main sys¬ 
tem we focus on the AGNR systems directly relevant for 
the lifting experiments.!^ The results for the ZGNR are 
provided mainly as comparison and to look into the role 
of chirality and in particular effects rooted in spin polar¬ 
ization, and thus we now discuss these separately. 


A. Pristine armchair nanoribbons 

As representative of the AGNR class we have investi¬ 
gated a pristine AGNR with a width of VF = 7 dimers 
(7-AGNR) corresponding to a G-G edge distance of 7.5 
A (see Eig.[^. It presents a direct semi-conducting band 
gap Eg due to the lateral confinement and can be classi¬ 
fied as a “large-gap ribbons” since p = 2 is an integer in 
the relation W = 3p 1.!^ We obtain Eg ^ 1.3 eV at the 
present level of approximation (DET-GGA and SZP ba¬ 
sis set), as seen from the electronic band structure shown 
in Fig.j^b). This value is smaller than those estimated 
experi mental ly {Eg « 2.3-2.6 eV for a flat AGNR on 
Au(lll)M^l^ and Eg ^ 2.7 eV for an AGNR suspended 
between surface and STM-tipP^ due to the underesti¬ 
mation of electron-electron interactiorP^ which plays an 
more important role in quasi one-dimensional GNRs com¬ 
pared to pristine graphene. Dielectric screening from the 
substrate also influences significantly the actual gap size: 
a band gap of 3.2 eV for a 7-AGNR was found to be 
lowered to 2.7 eV on a hexagonal boron-nitride (hBN) 
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(a) AGNR (b) ZGNR 
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FIG. 4: (Color online) lETS signals as a function of gate 
voltage for (a) pristine AGNR (4 vibrating unit cells) and 

(b) pristine ZGNR (6 vibrating unit cells). Vertical dashed 
lines are guides to the eye indicating the energy of the most 
contributing vibrational modes. Specific lETS signals for the 

(c) AGNR and (d) ZGNR at selected gate voltages marked 
with horizontal dashed lines in panels (a) and (b). Broad¬ 
ening originates from temperature T = 4.2 K and a lock-in 
modulation voltage Vms = 5 mV (except for the thin red lines 
in the lower panels with V-ms = 0 mV ). 

substrate using GW calculations S3 similar to the lower¬ 
ing calculated for a 7-AGNR on Au(lll).SSlln general we 
expect that underestimation of band gaps would mainly 
amount to a simple scaling the Fermi level position within 
the gap. 

We first discuss the effect of the finite size of the dy¬ 
namical region in our treatment. Figure [^a) shows how 
the lETS signals for the AGNR (at fixed gate voltage 
Vg = 0.0 V) vary as a function of the size of the dy¬ 
namical region, ranging from 1 to 6 unit cells. For easy 
comparison, the data are normalized by the number of 
vibrating unit cells. As the signal amplitudes in this rep¬ 
resentation are roughly constant we conclude that the 
absolute lETS simply scale linearly with the active e-ph 
coupling region. Consequently, the magnitudes in lETS 
may thus provide insight into the active scattering re¬ 
gion in actual experiments. Further, as we find that both 
lETS amplitude and shape is well converged with 4 vi¬ 
brating unit cells, we fix the dynamical region to this size 
in the following analysis. 

The computed lETS signals for the AGNR as a func¬ 
tion of varying gate voltage are shown in Fig. ga) as 
a density plot. Specific lETS spectra at selected gate 


FIG. 5: (Golor online) (a) Gomputed phonon band structure 
for the pristine, infinite AGNR (kx is in units of inverse unit 
cell length). The magnitude of the red, green and blue bands 
(corresponding to the three vertical lines in Fig.[^c)), is pro¬ 
portional to the signal size weighted overlap, {FnkiVc = OV) 
in Eq. between the repeated band vector and modes 

with frequencies huj > 180 meV, 180 > hw > 162 meV and 
hjj < 162 meV for red, green and blue, respectively. The 
red band is scaled by 0.2 compared to blue and green, (b-e) 
Selected phonon band modes at F for the infinite structure 
which, according to the projection, characterize the active 
lETS modes. Units of fc? 


voltages are shown in Fig. [^c) for both the intrinsic 
part (temperature broadening at T = 4.2 K) as well 
as that one would observe employing the experimental 
lock-in technique (additional broadening due to a mod¬ 
ulation voltage of V^rns = 5 mV). We find that for the 
AGNR there are generally two well-defined lETS signals 
appearing around 169 and 196 meV, corresponding to 
the D- (ring breathing) and G- {E 2 g phonon) modes, re¬ 
spectively, also observed in Raman spectroscopy.^^. The 
D-signal also has a shoulder with a local maximum at 
159 meV with contributions from several modes. These 
three distinct features are indicated with vertical lines 
in Fig. Qa,c). Shifting Ep inside the gap region with 
a relatively small gate voltage \Vg\ ^ 0.5 V does not 
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affect the lETS appreciably. However, when Ep comes 
close to the conduction band of the AGNR the signal in¬ 
creases by a factor of five and a small peak-dip feature 
appear similar to the on e repo rted for gated benzene- 
dithiol molecular contacts-l^^l^ Upon further gating into 
the conduction band the lETS signals undergo a sign re¬ 
versal (from peaks to dips) as the transmission increases 
beyond approximately 0.5 for the involved channels.!^ 
Similar effects are also found by gating into the valence 
band of the AGNR. 

We can easily identify the most important vibrational 
mode vectors vx for the lETS from the two amplitudes 
I7 aI and \kx\ given in Eqs. These modes can 

further be analyzed in terms of the phonons in the in¬ 
finite AGNR. To do so we introduce the measure F^k 
representing the overlap between modes in the finite dy¬ 
namical cell and the phonon band modes weighted by the 
size of the lETS signal, 


Fnk{VG)=E |7a( Vg)| ^Unk (^1, 


Ak 


AN-l)k^ 


■ Vx 


( 8 ) 

where Unk is the phonon band mode indexed by n, and Vx 
is the modes in a finite N primitive cell long dynamical 
region index by A. 

The projections EnkiVc = OV) are depicted as widths 
of the phonon bands in Eig. ^a), where the red, green 
and blue colors refer to modes with frequencies in the 
ranges hw > 180 meV, 180 > hw > 162 meV, and hw < 
162 meV, respectively. In total four bands contribute to 
the lETS signal corresponding to the four signals seen in 
the intrinsic part of the lETS spectrum in Eig. Sc). The 
corresponding T-point phonon modes inside the primitive 
cell for the infinite ribbon are shown in Eig. |^b-e). 


B. Pristine zigzag nanoribbon 

We next turn to our results for the pristine ZGNR 
shown in Eig. [^e). It has a width of IT = 4 zigzag 
“chains” (4-ZGNR) corresponding to a C-C edge distance 
of 7.26 A. The breaking of sublattice symmetry for the 
ZGNR and lack of pseudo-phase result in different selec¬ 
tion rules for the matrix elements and difference in for 
example Raman signals.^^The ZGNR generally presents 
spin-polarized edge states exhibiting a small band gap at 
the DET level,in our case ^ 0.6 eV (we note that this 
gap disappears in simpler tight-binding description^ or 
spin-degenerate DET calculations). The spin-polarized 
edge states play the major role for the conduction, see 
the spin-down eigenchannels visualized in Eig. [^e-h). 

Since the edge states break the mirror symmetry 
with respect to the middle of the ribbon, there are 
fewer symmetry-forbidden inelastic transitions between 
the scattering states for the ZGNR. Thus, we expect a 
wider range of modes to contribute to the lETS signal as 
compared to the AGNR case. Indeed this is in agreement 
with the findings shown in Eig.j^b) and Eig.|^b,d). The 
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FIG. 6: (Color online) (a) lETS signals for the pristine 

ZGNR (6 vibrating unit cells). The black lines correspond to 
spin-degenerate calculations while the red lines are the spin- 
up components of spin-polarized calculations. Broadening 
originates from temperature T = 4.2 K and modulation volt¬ 
age Urms = 5 mV (full lines) or Ums = 0 mV (dashed lines), 
(b) Electronic transmission from spin-degenerate calculations 
with varying electrode broadening describing the coupling to 
the metal contacts, 77 = 0,0.1,1 eV (see also Fig. dig) for the 
corresponding spin-polarized case). 



FIG. 7: (Color online) The phonon band structure of the 
ZGNR, {kx is in units of inverse unit cell length), together 
with the F-point modes. The widths of the red bands are 
proportional to the weight function E(0V) (Eq. ® ), while the 
widths of the blue bands are proportional to E(0V)+E(0.5V). 


greater number of modes contributing to the lETS for the 
ZGNR results in broader signals with similar magnitudes 
as compared to the lETS for AGNR. As for the AGNR 
case the lETS signal is well converged with a dynamical 
region consisting of 6 vibrating unit cells [Eig. [^b)]. 

Eor ZGNRs the ring breathing is forbidden by symme¬ 
try, thus the lETS is generally characterized by trans¬ 
verse and longitudinal modes. To explore the impact 
of spin-polarization on the ZGNR-IETS we compare in 
Fig.i the results from both spin-degenerate and spin- 
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(a) Clean 


(b) IH-edge 


(c) 2H-edge 








FIG. 8: (Color online) Top and side views of the dynamical region describing the various AGNR defect structures, (a) Pristine 
AGNR. (b) One extra H atom on one of the edges, (c) Two extra H atoms on one of the edges, (d) One H atom replaced by 
a F atom, (e) Dehydrogenated edge where 4 H atoms have been removed from each side, (f) One broken C-C bond, (g) Two 
broken C-C bonds, (h) Four broken C-C bonds, (i) Cu adatom in a hollow site on the edge. 


polarized calculations. Without gate voltage (Vg = 0 
V) the lETS display opposite signs due to the spin- 
induced gap. Only a single peak contributes to the 
spin-degenerate lETS while several peaks contribute to 
the spin-polarized lETS. Even if the ZGNR is tuned 
by Vg = 0.5V to become metallic and the two treat¬ 
ments then show the same overall sign in lETS, the spin- 
polarized lETS persists to show a much richer structure. 
This difference suggests that lETS could be a way to 
indirectly observe spin-polarized edge states. 

Projecting the modes contributing to the lETS onto 
the phonon band modes further underlines how several 
bands with different symmetries contribute to the spin- 
polarized lETS, while only a couple of bands contributes 
to the spin-degenerate lETS, see Eig. Again we use 
Eq. @ for this characterization, where the overlap for 
Vg = 0.0 V corresponds to the red color and the overlap 
for Vg = 0.5 V corresponds to the difference between the 
blue and red color in Eig. respectively. It is clear that 
spin-polarization permits more modes to contribute to 
the lETS. In contrast to the spin-degenerate case, where 
the symmetric electronic states (with respect to the mid¬ 
dle of the ribbon) only can couple to the symmetric vi¬ 
bration modes, the symmetry lowering of the electronic 
states by spin-polarization opens up also for scattering 
also via odd modes. 


IV. DEFECTED GRAPHENE NANORIBBONS 

In this section we address the modification and new 
signals in lETS that arise due to various defects in the 
GNR. Regardless of the fabrication method, defects will 


inevitable occur. For example, if the AGNRs are syn¬ 
thesized from a precursor molecule, involving heating 
and dehydrogenation, as reported by Cai et al\^ and 
Blankenburg et a/. there is a chance that the reaction 
is incomplete and some of the C-C bonds between the 
precursor molecules do not form. Also there is a chance 
that a part of the final AGNR will have dehydrogenated 
edges or are passivated by two hydrogen atoms. Finally, 
defects may be introduced on purpose by locally dosing 
a high current from the tip of a STM.l^ 


A. Defects in AGNRs 

In Fig. we show the structures of pristine AGNR 
along with 8 different defect configurations which we have 
considered. These include four defects in the edge passi¬ 
vation as follows: A single edge side with an extra hydro¬ 
gen atom [iH-edge, Fig. ^b)], two edge sides with each 
an extra hydrogen atom [2H-edge, Fig. |^c)], one hydro¬ 
gen replaced by a fluorine atom [iF-edge, Fig.j^d)], and 
a dehydrogenated edge with 4 hydrogen atoms removed 
from each side [8H-free, Fig.j^e)]. We have also consid¬ 
ered defects in the atomic structure in the form of one, 
two, or four broken C-C bonds [iC-broken, 2C-broken, 
4C-broken, Fig. |^f)-(h)] as well as a Cu adatom on the 
AGNR [Cu-adatom, Fig. [^i)]. For all these systems the 
entire dynamical region was relaxed, i.e., the parts of the 
AGNRs shown in Fig. 

Defects may influence the lETS signal in two ways. 
First, a defect can have a direct impact by changing the 
vibrational degrees of freedom. In order for the change in 
the vibrational spectrum to give a signal in the lETS, the 
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FIG. 9: (Color online) Electronic properties of the AGNR 
structures shown in Fig. The total transmission is shown 
with black lines. The ratio T/Ti, where Ti is the transmission 
originating from the most transmitting eigenchannel is shown 
with green dashed lines (this ratio gives a lower bound to the 
number of contributing eigenchannels). The DOS for the C 
atoms in the dynamical region is shown with red lines (offset 
by 3 units). 


new vibrations must couple to the current, and preferably 
have frequencies which do not coincide with ones already 
giving lETS signals for the pristine ribbons. Second, a 
defect can substantially change the electronic structure 
and thereby have an impact on the e-ph couplings as¬ 
sociated with the active modes or even the transmission 
eigenchannels of the pristine ribbons, e.g., changing a 
peak in the lETS to a dip (and vice versa) or enhancing 
asymmetric contributions via Eq. 

The electronic properties of the pristine AGNR is 
shown in Fig.j^a). The carbon DOS projected to the de¬ 
vice region (red curve) reveals a gap as expected from the 
band structure [Fig. [^b)], which is significantly broad¬ 
ened from the coupling to the metallic electrodes. The 
two valence and two conduction bands in the considered 
energy range naturally explain that the total transmis¬ 
sion (black curve) is bound below a value of 2. Further, 
the ratio T/Ti < 2 (green dashed line), measuring the 
minimum number of contributing channels where T\ is 
the transmission of the most transmitting eigenchannel, 
shows that both channels play a role for the transport, 
at least away from the edges of the direct band gap. 
Measurements of shot noise may provide insights into 


this effective number of conductance eigenchannels. 

We can now discuss how the different defects modify the 
electronic properties. From Fig. [^b)-(i) we notice that 
not all defects change the elastic transmission, and fur¬ 
thermore, a change in elastic transmission needs not be 
unique for a specific defect. 

Instead, lETS may provide a additional fingerprint in 
the current that can be used to identify the type of de¬ 
fect. Figure pT] shows the computed lETS as a function 
of gate voltage for the 8 different defects. As for the 
clean structure, the two peaks at 169 and 196 meV corre¬ 
sponding to the D- and G- Raman modes are dominant 
for a range of gate values for all the structures. Another 
feature, which is present in all the systems, is the appear¬ 
ance of several signals close to the band onsets. In the 
following subsections we discuss in more detail the trans¬ 
port characteristics with the different types of defects in 
AGNRs. 


1. Edge passivation 

Considering defects in the edge passivation [Fig. [^b- 
e)] the gap in the transmission is essentially unchanged 
[Fig. [^b-e)], except for the IH-edge structure where a 
zero-energy resonance appears in the DOS and trans¬ 
mission [Fig. [^b)]. This new peak can be attributed 
to tunneling via a mid-gap state which appears due to 
the local breaking of sub-lattice symmetry.^ Thus, if a 
H atom is added to the neighboring C atom [2H-edge, 
Fig. I^c)] the peak disappears [Fig. |^c)]. The addition 
of one or two H atoms on the same side also results in the 
closing of one transmission channel between the valence 
and conduction bands as shown in Fig. [^b,c). Concern¬ 
ing the vibrational degrees of freedom, the addition of 
extra hydrogen to the edge results in new vibrational 
modes around 330 meV for IH-edge and around 343 and 
353 meV for 2H-edge, clearly outside the bulk phonon 
band (ranging up to ^ 200 meV) of pristine AGNR.l^ 

Comparing the lETS in Fig. [T^a-c) we find that only 
IH-edge gives a signal which differs significantly from 
the pristine case. Figure [^k) shows specific lETS for 
selected gate voltages for IH-edge. Here, at Vq = 0.2V 
(top green curve) we see how new signals appear at large 
voltages: For positive bias polarity two signals appear at 
330 and 365 meV, respectively, while for negative bias 
polarity only an asymmetric signal around —365 meV is 
present. The signal at 330 meV is due to vibrations of the 
H2 [Fig.[^b)], while the signal at 365 meV [Fig.pTj^a)] is 
due to the H atom on the neighboring C atom. Further, 
the amplitude of the signals around 169 and 196 meV is 
also found to depend on bias polarity. 

Gating onto the zero-energy resonance for IH-edge the 
lETS signal [middle red curve in Fig. m^)] is dom¬ 
inated by large asymmetric signals for low energy vi¬ 
brations due to the contribution from kx and Eq. Q. 
We note that kx changes sign with bias polarity for this 
approximately left-right symmetric structure. This can 
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FIG. 10: (Color online) (a-i) lETS as a function of gate voltage Vg for the pristine and defected AGNR structures shown in 
Fig. [^. (j-o) lETS for six selected structures at three specific gate values (dashed horizontal lines in panels a-i). The curves 
are offset with the most negative gate value at the bottom (black curves) and the most positive at the top (green curves), (j) 
Clean AGNR at gate values Vg = —0.3, 0.0, and 0.8 V. (k) IH-edge at Vg = —0.3, 0.0, and 0.2 V. (1) 8H-free at Vg = —0.3, 
0.0, and 0.6 V. (m) IC-broken at Vg = —0.3, 0.0, and 0.3 V. (n) 4C-broken at Vg = —0.3, 0,0, and 0.3 V. (o) Cu-adatom at 
Vg = —0.3, 0.0, and 0.3 V. Dotted vertical lines are guides to the eye of characteristic lETS signals corresponding to the modes 
in Fig. 


be seen from the red lETS curve in Fig. 10 ^k) which is 
roughly an odd function of the bias voltage. In close 
proximity of the zero-energy resonance a characteristic 
“X-shape” is observed in the gate-dependent lETS, while 
away from it the signals approach that of the pristine 
AGNR [Fig.j^b)]. 


Substituting a H atom with a F atom (IF-edge) is 
seen to have virtually no effect in the lETS of Fig.[^d). 
This suggests that a significant change in the chemical 
composition directly involving the 7r-electronic system is 
required in order to obtain a signal although the vibra¬ 
tions are influenced by the heavier passivation. 


Such a significant change in the passivation occurs 
for instance by removing four H atoms on each side 
(8H-free), giving rise to four very narrow peaks in the 
DOS around the conduction band, [Fig. [^e)]. These 
correspond to very localized dangling-bond states on 
the dehydrogenated dimers and therefore do not show 
up in the transmission. However, the dehydrogenated 
edges give rise to localized vibrations outside the range 
of the pristine vibrational spectrum.!^ The in-phase vi¬ 
bration of the dehydrogenated C dimers at the armchair 
edges [Fig. &)] gives rise to an extra lETS peak at 


244 meV [Fig. pT]^!)] matching the H-free mode measured 
by Raman .1^ We find that this signal is robust as it ap¬ 
pears in the whole range of gate values. When gating 
into to the valence band a new signal appears around 
43 meV [Vg ~ —0.8 V in Fig. p!Q[ e)] originating from a 
low energy edge vibration [Fig. |ll[g)]. 


2. Structural defects 


The electronic transmission in GNRs is mediated by 
the carbon tt system. Thus if a C-C bond fails to be 
formed during GNR synthesis or if it is broken again 
at a later stage, a large effect can be expected for the 
electronic conduction properties. This impact is indeed 
revealed in Fig. [^f-h). Breaking one or two bonds re¬ 
sults in the formation of two in-gap states which, broad¬ 
ened by the electrodes, make the gap appear smaller. 
The lETS signals for the IC-broken and 2C-broken in 
Fig. 10^f,g,m) have the same two signals at 169 and 196 
meV as for the clean ribbon. However, the relative am¬ 
plitudes are interchanged such that the "D”-peak is now 
slightly more intense than the ”G"-peak. 
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FIG. 11: (Color online) Visualization of the most contributing defect-induced vibrational modes to the lETS signals indicated 
by vertical lines in Fig. [To|j-o). (a-b) The two hydrogen signals for IH-edge. (c) Localized edge mode at the carbon dimers for 
the 8H-free. (d) Delocalized edge mode for the 8H-free. (e) Hydrogen mode from the zigzag edge of 4C-broken. (f-g) Defect 
modes for 4C-broken. 


Breaking four C-C bonds [4C-broken, Fig. |^h)], re¬ 
sulting in constrictions of single C-C bonds, totally alter 
the DOS which is now dominated by three sharp peaks 
as seen in Fig. [^h). The corresponding lETS signals 
are shown in Fig. IT^h,n). In the proximity of the zero- 
energy resonance a broad range of signals at low vibra¬ 
tional energies appears (red curve in panel n) as well as 
a characteristic “X-shape” in the gate plot (panel h) sim¬ 
ilar to that of IH-edge. Gating away from the resonance 
we observe two additional robust lETS signals at 27 and 
50 meV resulting from vibrations localized at the defect 
[Fig.[TiId,e)|. 


3. Adatoms 

Transition metals are typically used for growth of 
graphene or as a substrate for the bottom-up synthe¬ 
sis of GNRs. Thus it is of interest to consider the effect 
of adatoms of this type on GNRs. A Cu adatom on 
graphene adsorbs preferentially in the on-top position.l^ 
However, positioning Cu such that it breaks the axial 
symmetry of our AGNR, we find that it is most stable in 
a hollow site at the edge JCu-adatom, Fig.j^i)]. The DOS 
and transmission in Fig.l^i) reveal a n-type doping effect 
shifting Ep close to the conduction band while leaving 
the two transmission channels inside the gap relatively 
intact. 

For the pristine GNR the e-ph couplings of the out-of- 
plane vibrations are suppressed due to the symmetry of 


the TT-orbitals. However, around the onset of the conduc¬ 
tion band the lETS signals in Fig. 10p,o) is dominated 
by large asymmetric signals with significant contributions 
from out-of-plane phonons. These modes come into play 
due to breaking of the planar symmetry by the adatom. 
Also note that by gating of Ep within the gap these sig¬ 
natures of the adatom disappear, cf. the lower black curve 
in Fig.lMo). 


B. Defects in ZGNRs 

Let us next consider a series of defects for the zigzag 
graphene nanoribbon. In Fig. we show the atomic 
structures of pristine ZGNR along with 8 different de¬ 
fect configurations. We consider the following defects in 
the edge-passivation: A single edge with an extra hy¬ 
drogen [IH-edge, Fig. [T^b)], one hydrogen is replaced 
by either a F atom [iF-edge, Fig, p^c)], an OH group 
[lOH-edge, Fig. [^f)], or a NO 2 group [lN02-edge, 
Fig. [^g)]- We also consider defects in the form of 
a Cu adatom [Cu-adatom, Fig. [T^d)] or a Li adatom 
[Li-adatom. Fig. [T^e)]. Finally, we also study the ef¬ 
fect of a structural defect in form of a 57 reconstruc¬ 
tion [R57, Fig. p!^h)] and a substitutional defect where 
a C atom next to the edge is replaced by a Si atom 
[Si-substitute. Fig. [T^i)]. For all these systems the 
entire dynamical region was relaxed, i.e., the parts of 
the ZGNRs shown in Fig. [fusing spin-polarized treat¬ 
ments. The spin degrees of freedom a generalizes 
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FIG. 12: (Color online) Top and side views of the dynamical region describing the various ZGNR defect structures, (a) 
Pristine ZGNR. (b) One extra H atom on one of the edges, (c) One H atom replaced by a F atom, (d) Cu adatom in a hollow 
site on the edge, (e) Li adatom in a hollow site on the edge, (f) One H replaced by a OH group, (g) One H replaced by a NO 2 
group, (h) Structural defect (R57). (i) Substitutional Si defect next to the edge. 
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FIG. 13: (Color online) Electronic properties of the ZGNR structures shown in Fig. with the spin-up/down components 
in the left/right panel. The spin-resolved total transmission is shown with black lines while spin-averaged total transmission 
is shown with thin blue lines. The ratio where T[ is the transmission originating from the most transmitting spin 

eigenchannel, is shown with green dashed lines (this ratio gives a lower bound to the number of contributing eigenchannels with 
spin a). The spin-resolved DOS for the C atoms in the dynamical region is shown with red lines (offset by 3 units). 
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7a and [Eqs.0-@] corresponding to two indepen¬ 
dent spin channels, which in general can have quite dif¬ 
ferent amplitudes and even opposite sign. The observable 
lETS would simply be the sum of these two components 

{dm+dm)/{dvh+dvh). 

Similar to the AGNR case, the electronic properties in 
the device region with the different impurity configura¬ 
tions for the ZGNR, now spin resolved, are summarized 
in Fig. The lETS of pristine ZGNR was already dis¬ 
cussed in Sec. |IIIB] and below we continue describing the 
lETS fingerprints for the various defects. 


1. Edge passivation 

As commented above, the broader lETS signals of pris¬ 
tine ZGNR [Figs. [3][^ (as compared with AGNR) can 
be understood from the breaking of the axial mirror 
symmetry and hence fewer symmetry-forbidden inelas¬ 
tic transitions. These broader signals may in general 
make the detection of defect signatures more difficult. 
For IH-edge [Fig. P^b) ] the lETS resembles that of the 
pristine ZGNR [Fig. |13[a)] inside the gap. However, gat¬ 
ing into the valance band [black curve in Fig. p^k)] the 
edge states start to extend into the middle of the ribbon, 
partially restoring mirror symmetry, and thus resulting 
in part of the pristine ZGNR signals to disappear. 

Here an extra signal appear due to edge-modes in the 
frequency range 194 to 199 meV with the most contribut¬ 
ing mode at 196 meV as shown in [Fig. [T^a)]. The re¬ 
sulting lETS signal can clearly be seen in the bottom 
curve in Fig. p^k). As for the AGNR substituting a hy¬ 
drogen with a fluorine atom has a very limited effect on 
the electronic properties and the lETS signal. 

Substituting a hydrogen with an OH group, accord¬ 
ing to Fig. p!^f) and (o), have only a small effect on 
the spin down electrons, while it shrinks and add ad¬ 
ditional structure to the gap for the spin up electrons. 
For the spin up electrons there is a small peak inside the 
gap which gives rise to a large asymmetric lETS signal 
around Vq = —0.2V in Fig, p^o) lower curve, compared 
to the pristine case. The most contributing mode to the 
asymmetric lETS signal is shown in Fig.[T^b). However, 
there is no clear signature of the OH group itself. In the 
same manner the substitution with a NO 2 group removes 
the gap in the electronic properties without leaving any 
direct fingerprint of the NO 2 group in the lETS signal. 


2. Adatoms 


contributions. For a gate value of Vg = —0.2 V, the lETS 
signal is dominated by spin down electrons. Due to the 
finite width of the in-gap peak, in the spin down trans¬ 
mission, the low frequency phonons {hw <0.1 meV) give 
rise to back scattering while the high frequency phonons 
{huj > 0.1 meV) result in forward scattering. Thus, the 
low and high energy signals have different signs as can be 
seen from Fig. [T^l). Interestingly, the low energy signal 
primarily consists of symmetric contributions from out- 
of-plane modes [Fig. p!^c)]. Replacing the Cu adatom 
with Li, the transmission and DOS, shown in Fig. 13 ^e,n), 
reveals a spin dependent n-type doping effect, where Ep 
is shifted the most for spin down. However, no in-gap 
peak is seen as for Cu and the lETS show no clear sig¬ 
nature of the Li atom. 


3. Structural defect 

The formation of a R57 reconstruction results in peaks 
in the DOS in the device region, just above Ep for spin up 
[Fig.p^h)] and just below Ep for spin down [Fig.[l^q)]. 
The R57 breaks the symmetry both in the vibrational 
and electronic structure allowing for lETS signals from 
a wider range of vibrations, resulting in broader peaks, 
as seen from Fig. Elh) and Fig. [l^n). One of the con¬ 
tributing modes is localized at the border between the 
pentagon-ring and middle of the ribbon at hw = 204meV 
[Fig. p!^d)]. This localized mode yield a relatively small 
signal compared to the other signals, however, contrary 
to the other modes the localized mode is not expected to 
be broadened if the coupling to phonons away from the 
dynamical region is taken into account. The breaking 
of symmetry in the electronic structure also give rise to 
difference signals for the two bias polarities. 


4 . Substitutional impurity 

Substituting a carbon with a silicon atom leads to an 
out-of-plane buckling, see Fig. [T^i). However, both sil¬ 
icon and carbon have an electronic structure, and 
the electron transmission is basically similar to the pris¬ 
tine. On the other hand, the buckling give rise to low en¬ 
ergy peaks in the lETS signal originating from the e-ph 
coupling to the out-of-plane modes [Fig. [T^e)]. Gating 
close to the band edge of the conduction band gives rise 
to different sign of the signals at low and high vibrational 
energies, as seen from the top curve in Fig. [l^o). 


As for the AGNR we consider the effect of adatoms. 
For the Cu adatom the transport gap shrinks for the spin 
up electrons while there is an in-gap peak for the spin 
down electrons, cf. Fig.p^m). Thus, for some gate values 
the lETS signals reflect that the spin down electrons will 
back scatter while the spin up electrons will be forward 
scattered, and the observed signal is then the sum of these 


V. CONCLUSIONS 

In summary, we have investigated lETS signals in 
symmetrically contacted armchair and zigzag graphene 
nanoribbons, considering both pristine as well as a selec¬ 
tion of defected configurations under varying charge car¬ 
rier conditions. For the clean AGNR inelastic tunneling 
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FIG. 14: (Color online) (a-i) Spin-averaged lETS as a function of gate voltage Vg for the pristine and defected ZGNR 
structures shown in Fig.[^ (j-o) lETS for six selected structures at three specific gate values (dashed horizontal lines in panels 
a-i). The curves are offset with the most negative gate value at the bottom (black curves) and the most positive at the top 
(green curves), (j) Clean ZGNR for gate values Vg = 0.0, 0.4, and 0.8 V. (k) IH-edge at Vg = —0.6, 0.0, and 0.4 V. (1) 
Cu-adatom at Vg = —0.2, 0.0, and 0.8 V. (m) lOH-edge at Vg — —0.2, 0.0, and 0.8 V. (n) R57 at Vg = —0.8, 0.0, and 0.2 V. 
(o) Si-substitute at Vg = —0.3, 0.0, and 0.3 V. Dotted vertical lines are guides to the eye of characteristic lETS signals 
corresponding to the modes shown in Eig. |15| 


gives rise to two distinct peaks in the lETS spectrum at 
169 mV and 196 mV corresponding to the D-and G-modes 
of Raman spectroscopy, respectively. By connecting the 
lETS signals to the phonon band structure, we have clar¬ 
ified how only a single band contributes to the ”G-mode” 
while three bands contribute to the broader "D-mode". 
Concerning defects in AGNRs we have shown how some 
leave lETS unchanged while others give clear signals. For 
instance, adding an extra hydrogen atom to a single edge 
side gives a clear signal for some gate values. This signal 
can be removed by adding another hydrogen atom to the 
neighboring edge side because the sub-lattice symmetry 
is restored. Further, exchanging a single hydrogen atom 
with a fluorine atom in the passivation does not result 
in any change in both the elastic and inelastic tunneling. 
However removing 8 hydrogen atoms leaving part of the 
edge on each side without passivation, gives a clear ro¬ 
bust signal throughout the investigated gate values. The 
signal, due to the vibration of the carbon dimers at the 
edge, has an energy around 245 meV making it easy to 
detect since it is outside the vibrational spectrum of the 
pristine ribbon. Breaking of one or two C-C bonds turns 
out to interchange the relative intensity of the "G"- and 
"D"-peaks. Breaking 4 C-C bonds gives rise to signals 


caused by the defect tilted out of plane. Lifting the sym¬ 
metry of the TT-electrons by adding a Cu-adatom allows 
the out-of-plane modes to contribute. 

For the ZGNR we find relatively broader lETS signals 
especially in the absense of a large gate voltage (Vg ^ 0 
V). Importantly, this is a consequence of the breaking 
of the axial mirror symmetry in the ribbon due to the 
presence of spin-polarized edges. Thus, by comparing to 
spin-degenerate calculations, we suggest that lETS can 
give an indirect proof of spin-polarization in zigzag rib¬ 
bons. On the other hand, the broader lETS features may 
make it difficult to identify the different defect signals re¬ 
ported in this paper. 

The presence of a R57-reconstruction also broadens the 
lETS by breaking both the electronic and vibrational 
symmetry. Substituting a carbon atom with a silicon 
atom makes the ribbon buckle, breaking the planar sym¬ 
metry, allowing the out-of-plane modes to contribute to 
the lETS. This suggests that lETS in principle could be 
used to gain information of the curvature of GNRs and 
other graphene-based structures. 

Finally, as an outlook we note that here we presented 
calculations on long, symmetrically contacted systems 
where there is a significant overlap with both metallic 
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FIG. 

Fig. 


15: (Color online) Visualization of the most contributing modes to the lETS signals indicated by vertical lines in 

j-o). (a) Edge mode at the edge with the extra hydrogen in IH-edge. (b) Mode contributing to the asymmetric signal 
(m) for the lOH-edge. (c) Out-of-plane mode for Cu-adatom. (d) Localized mode for R57. (e) Out o plane mode for 
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in Eig. 
Si-substitute. 


electrodes. It would be interesting to extend such a study 
also to the asymmetric situation where a point tunnel 


contact is made to one end resembling, say, the coupling 
a STM tipP 
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